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We have conducted two quantitative tests of predictions based on the Halperin-Lubensky-Ma 
(HLM) theory of fluctuation-induced first-order phase transitions. First, we explore the effect of 
an external magnetic field on the nematic-smectic-A (NA) transition in a liquid crystal. Second, 
we examine the dependence of the first-order discontinuity as a function of mixture concentration 
in pure 8CB and three 8CB-10CB mixtures. We find the first quantitative evidence for deviations 
from the HLM theory. 



One of the most important advances made in our un- 
derstanding of continuous phase transitions has been the 
modification of critical exponents, due to thermal fluc- 
tuations, from the values predicted by mean-field the- 
ory But thermal fluctuations have another effect, 
one that is less-well-understood theoretically and stud- 
ied only to a limited extent experimentally: when two 
order parameters (or an order parameter and a gauge 
field) are simultaneously present and interact with each 
other, the fluctuations of one may drive the phase tran- 
sition of the other first order. In high-energy physics, 
for example, an analogous situation occurs in the Higgs 
mechanism Q. In condensed-matter physics, over two 
decades ago, Halperin, Lubensky and Ma (HLM) ||] 
predicted that this could occur in two settings: the 
normal-superconducting phase transition in type-1 su- 
perconductors and the nematic-smectic-A (NA) tran- 
sition in liquid crystals. At the NA transition |3] [|, 
the first fluctuating field is the nematic order parame- 
ter Qij = S (SriiUj — Sij) /2, while the other order pa- 
rameter is the two-component smectic-A order param- 
eter ip. In principle, a complete theory must account 
for thermal fluctuations in the nematic order parameter 
magnitude S, the nematic director n, and the smectic-A 
order parameter ip. The de Gennes-McMillan mecha- 
nism [|| takes into account S fluctuations but is mean 
field in n and ip. The HLM theory, in addition, considers 
n fluctuations but is mean field in ip. Deviations from 
the de Gennes-McMillan mechanism have been quantita- 
tively established by Anisimov et al. Q and Cladis et al. 
H . Although their conclusions are consistent with the 
HLM mechanism, the experiments are at their resolution 
limits and the full implications of the HLM theory remain 
to be tested. 

In this paper, we conduct the first two quantitative 
tests of the HLM theory, both leading to significant dis- 
crepancies from HLM. (1) Mukhopadhyay et al. |J have 
made a detailed prediction for the external-field depen- 
dence of the HLM effect; we look for this field effect ex- 
perimentally in the cyanobiphenyl liquid crystal 8CB, as 
well as in mixtures of 8CB and 10CB. (2) The first-order 
discontinuity implied by the HLM theory (which can be 
characterized by the number to = (Tna — T*)/T*, where 
Tna is the transition temperature and T* is the spin- 
odal point) should also be sensitive to the mixture con- 



centration x in the 8CB-10CB system. Here, we make a 
quantitative study of to vs. x. 

In order to achieve the resolution required for these 
measurements, we designed a high-resolution real-space 
optical technique ("Intensity Fluctuation Microscopy") 
which is at least an order of magnitude more sensitive 
than previous techniques. The experimental set-up is 
briefly described next (see Ref. jl(|[ll| for more details). 

The liquid crystal is sandwiched between two glass 
cubes, both treated for unidirectional planar anchoring. 
The spacing between the cubes (typically d — 30 //m) 
is adjusted interferometrically, and thickness gradients 
across the entire cell are controlled to be 0.5 /im. The 
cubes are placed in a gradient hot-stage specifically de- 
signed for this experiment to fit an in-house transverse- 
bore electromagnet. The top and bottom of the sample 
cell are independently temperature controlled at temper- 
atures that straddle the NA transition temperature. 

We use a crossed-polarized digital CCD microscopy 
set-up to measure intensity fluctuations as a finely re- 
solved function of temperature. These intensity fluctu- 
ations are caused by director fluctuations which are a 
soft mode in the nematic phase but are suppressed in the 
smectic-A phase. Subtracting two successive, uncorre- 
lated "instantaneous" images (a 1.8 mm x 1.4 mm area of 
the sample being imaged onto a 1317 x 1024 pixel, 12-bit 
digital CCD camera) gives us a difference image 



A(x,y) = h{x,y) - I 2 (x,y), 



(1) 



which reflects static intensity fluctuations. The his- 
togram of the intensity differences in the nematic and 
smectic-A phase is well fit by a Gaussian probability dis- 
tribution. The primary contribution to the noise is pho- 
toelectronic shot-noise, which also has a Gaussian prob- 
ability distribution. The normalized variance, £ has the 
shot-noise variance subtracted and is normalized with re- 
spect to the transmitted light intensity: 



c-c 



((I) - Ioffset) 2 ' 



(2) 



where ( is the measured raw variance, (sn is the shot- 
noise variance calibrated in a blank field, ( ) refers to a 
spatial average over pixels, (I) is the average intensity 
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over both images and all pixels, and I ffset is an built-in 
intensity offset in the CCD camera. 

The variance £ decreases as the NA transition temper- 
ature, Tna, is approached from above. In the smectic-A 
phase, it is much smaller; in fact, £ is not measurably 
higher than the zero (shot-noise) level. Since the applied 
temperature gradient straddles the NA transition, the re- 
sulting difference image contains a hot nematic part and 
a cooler smectic-A part separated by a flat boundary. We 
divide our image into strips parallel to the NA boundary, 
which is typically flat over the field-of-view. By calibra- 
tion, we translate the spatial dependence of £ into a tem- 
perature dependence. Each image that contains an NA 
interface is self- calibrating in that over small times, the 
gradient remains constant while the overall temperature 
fluctuates slowly. The r.m.s. fluctuations in temperature 
over long times (« I hour) was roughly 0.15 mK and was 
less than 0.05 mK over the duration of a measurement. 

The experimental technique described above was first 
used to study the critical behavior close to Tna- Pure 
8CB is a material that exhibits a NA transition that 
is indistinguishable from second order when studied by 
traditional high-resolution calorimetry. Dynamical mea- 
surements by Cladis et al. || have suggested that this 
transition is first order. We confirm here (see also |10|) 
that this transition is unambiguously first order. 
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FIG. 1. Pure 8CB, d = 30^m, 0.59 ± 0.02K/cm. Here, 
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5.0±0.5 mK. 



1.6 ±0.1 x 10" 



Results are shown with 



and without an external applied field. 

The fluctuation- vs. -temperature profile for 8CB is plot- 
ted in Fig. |l| (black data curve). There are three distinct 
sections to the graph: 

1. Smectic: Here, C = 3.1±3.0 x 10 -6 , and the smectic 



fluctuation level is flat and indistinguishable from 
the shot-noise background. 

2. Nematic: The data here fit a power law of the form 

w Q + Wl (T - w 2 ) W3 (3) 

with w = 3.1±3.0x 10~ 6 , w x = 1.43±0.04 x 10~ 5 , 
w 2 = 5.0 ± 0.5 mK, and w 3 = 0.50 ± 0.05. The fit 
(shown in Fig. 0) was done over the largest tem- 
perature range that still kept the residuals (top 
curve in Fig. |]) flat, wq was held to the smec- 
tic level, while w% is the exponent of the power 
law. W2, the divergence point of the power law, is 
identified with the spinodal temperature T* . How- 
ever, it is a phantom divergence point, as the actual 
phase transition is located at the interface, i.e., at 
T-T NA = 0mK @ 

3. Interface: Ideally, the interface on Fig. [j] would be a 
vertical line. However, a temperature gradient nor- 
mal to the sample plane smears out the interface. 
Varying temperature moves the position of the in- 
terface in the field-of-view, but the magnitude of 
the smearing is stable to within the statistical er- 
rors in the data. In practice, we fit the entire data 
to a convolution of the power law form in Eq. ^ 
and a linear temperature gradient. This fit differs 
from the simple power law fit (also shown) only in 
the interfacial region. The position of the interface 
is unambiguously determined by the fit and does 
NOT coincide with the zero point. The top curve 
in Fig. |l| shows the residuals for such a fit. 

We tested for the effect of sample thickness d on the 
measurements. We find no systematic dependence on 
sample thickness for d > 7.5 /im. For smaller d, we find 
a marked increase in to, which we ascribe to the effect of 
a small residual rotation (« 0.2°) between the anchoring 
direction of the top and bottom plates. The twist created 
by this rotation is strongest for small d. The implications 
of this will be examined in a longer article []l3|| . 

Having established the experimental technique, we first 
study the critical effects of suppressing director fluctua- 
tions with an external magnetic field applied along the 
easy axis of the nematic. Magnetic-field effects in the 
nematic phase have been studied in the last two decades 
0,^5|. Scale fields to change Landau parameters at the 
NA transition are large; the mean-field effects of such 
large fields have been studied by Lelidis and Durand jl6| . 
Recently, Mukhopadhyay et al. Q have shown that based 
on the HLM mechanism, one expects modest fields to 
suppress the director fluctuations that drive a mean-field 
second-order transition first order. They predict a non- 
analytic dependence of to on the field H and a critical 
field, H c , above which the transition is second order. In 
particular, they predict 



H r = Hn t. 



o t 



1/2 



(4) 
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where Hq is a scale field that depends on known mate- 
rial parameters. For 8CB, H as 3500 T and from Fig. ||, 
t = 1.6 x 10~ 5 : this implies a critical field of H c * 10T. 

The magnetic field effect at the transition was probed 
with a field that could be varied from 0T to 1.5 T. Con- 
trary to initial expectations, no suppression of the fluc- 
tuations was seen within 50 mK of the transition point. 
Two data curves are shown in Fig. [|. The data curves 
(black representing zero field, and grey representing mea- 
surements at a field of 1.2 T) practically fall on top of each 
other. Likewise, we see no shift in t . 

What do the magnetic field studies on a 0-1.5T field 
tell us about the critical field of ss 10 T predicted via the 
HLM scenario? A weaker effect implies a larger critical 
field. 

To establish an lower bound for H c , we look at the 
small- 1 H | limit of the HLM prediction, 

t (H)/t nl-mH/H c , (5) 

where m = (2.17 ± 0.01) is the slope of the linear sup- 
pression predicted in Ref. M. In our case, our confidence 
level for to is roughly 10%; i.e., the largest suppression 
of to that we would not observe is to{H)/to = 0.9. This 
puts a lower bound on the value of H c . The maximum 
field used is H = H max = 1.5 T. Thus, for 8CB, we find 



c " l-t (H)/t 
_ (2.17 ±0.01) 1.5 T 
~~ 09 

«33T. (6) 

Thus, the lower bound for the critical field in our mea- 
surements is three times the predicted value. 

We also measured the effect of a magnetic field in two 
mixtures, 0.18 and 0.41 mole fraction 10CB in 8CB. In 
neither of these samples do we observe a suppression of 
fluctuations close to Tma- 

One might worry about the validity of a null result. We 
do, however, see the expected non- critical suppression of 
the fluctuations. We find that this depression falls off 
quadratically with increasing field: 

C = Co-Sff# 2 - (?) 

This non-critical reduction in the field-effect on ap- 
proaching Tma reflects quite clearly the fact that the 
(twist and bend) elastic constants are getting stiffer as 
one approaches Tna- I n fact, since the magnetic coher- 
ence length is = ( y ) 2 jj , one reasonably expects that 
the field dependence is always a function of i.e., the 
coefficient must obey gn oc -4, where K is a diverging 
elastic constant. Fitting to the value of v = 0.5 that is 
consistent with other experiments and with our earlier fit 
in zero field, we find that the data clearly agree with this 
(see also J13[.) 

Our second experimental test of the HLM effect was 
the dependence of to on the concentration x of 10CB 



in 8CB, between x = and x = 0.44. This depen- 
dence has been studied previously by adiabatic calorime- 
try (Marynissen et al. E^]). Anisimov et al. deter- 
mined the Landau tricritical point (the LTP is the tri- 
critical point predicted by the de Gennes-McMillan the- 
ory) to be at x* — 0.41. The resolution in the adiabatic 
calorimetry is, however, limited to x > 0.3. jl8) To test 
the HLM theory unambiguously, we extended the range 
of mixture concentration to x — (pure 8CB). The non- 
zero latent heat for x < x* w 0.41 implies that their 
data are inconsistent with the de Gennes-McMillan the- 
ory but consistent with the HLM form, which contains 
a small negative cubic term (with a constant cubic coef- 
ficient B) that keeps the latent heat non-zero. In HLM, 
the effective mean-field Landau free energy is 

/ = \a\iP\ 2 - ii?H 3 + hj\^ + \e\^ (8) 

In Ref. [Q, the condition for coexistence (/ = and 
dF/dip = 0) was used to calculate the entropy change 
AS = L/T, in terms of the three Landau parameters A, 
B, and C. Close to the transition, A = at, where t is 
the reduced temperature, and a is a constant. The vari- 
ation of the quartic coefficient near the LTP is modeled 
by C = Co {x — x*). To make connection with our mea- 
surements of to, we calculate to = An a /a, where A n a, 
the quadratic coefficient at the transition, is a function 
of B, C, and E. 
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FIG. 2. (a) Fit of to data to the Anisimov parameters, with 
Co as a fit parameter. Top and bottom dashed curves show 
fits to to for v — 0.6 and v — 0.4, respectively. (The data 
points for these exponent values are not shown.) (b) Com- 
parison of latent heat data (Ref. 16) to to data (this work). 
Open circles: data taken from Ref. 16; filled circles: to data 
converted to equivalent latent heat. 
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The fit parameters from the analysis of Anisimov et al. 
are 



= -3/8(aC /E) = 0.993 
AS*/R = 0.0261 



(9) 



While the entropy change per mole AS/ R depends only 
on B, aCo and E, the discontinuity to depends on all 
four parameters separately: here, we use the two Anisi- 
mov parameters and fit Co and E. 

Each data point in Fig. ||(a) is obtained from the 
power-law fit described earlier. The exponent of the 
power law is fixed at v = 0.5. The error bars on a free 
fit are ss 0.1. In Fig. ||(a) we plot the data for a fit- 
ted exponent of v = 0.5 and the HLM fit. Also shown, 
for perspective, are HLM fits for to data with exponents 
v = 0.4 and v = 0.6. The values of Co in the three fits 
are -0.74 ±0.03, -0.79 ±0.03, and -0.88 ±0.03, respec- 
tively. (The best-fit value of E = 0.97 ± 0.03 for D = 0.5 
is used for all three fits.) In Fig. |2|(b], we compare the 
latent-heat data of Marynissen et al. Jlj] with that calcu- 
lated from our t data (using the Anisimov parameters). 
The discontinuity gets larger as we approach the Landau 
tricritical point. In particular, we find that the data near 
the LTP agree with the HLM prediction for reasonable 
values of the quartic coefficient Cq. However, for lower 
concentrations of 10CB in 8CB, there is a clear devia- 
tion from the HLM prediction, which underestimates the 
latent heat in this region of the phase diagram |1| . 

Two main results arise out of this study, and both 
constrain theories in the weakly type-1 region of the NA 
transition. The first is the absence of a magnetic field 
effect on the critical behavior in pure 8CB and in mix- 
tures. It was predicted 0] that an external field would 
suppress nematic director fluctuations and drive the tran- 
sition back to second order. Experimentally, we see a 
non-critical suppression of the fluctuations but observe 
no change in the critical behavior. In particular, the 
strength of the first-order transition, measured by the 
dimensionless number to, is unchanged. We thereby es- 
tablish a lower bound of 30 T on the critical magnetic 
field, three times higher than the predicted value. We 
conclude that the neglect of smectic fluctuations by HLM 
is not valid in this region of the phase diagram. 

Second, the variation of to in mixtures has yielded the 
surprising result that to on the weakly first-order side of 
the LTP is larger than the HLM predicted value, and 
not smaller as one might naively expect of a crossover 
to second-order XY behavior. This might be related to 
the observed local maximum in critical exponents seen 
at the 3DXY to tricritical crossover |Q. This result puts 
a strong constraint on any theory of the tricritical point. 

Both results in this work provide, for the first time, 
quantitative evidence for deviation from the implications 
of both the de Gennes-McMillan theory and the HLM 
theory. Our data show not only that the discontinuity 
is larger than that predicted by HLM theory but also 
that the field dependence is much weaker than predicted. 



These results provide a testing ground for any future the- 
ory that goes beyond HLM; the logical next step would 
be to account for the effect of smectic fluctuations. Fu- 
ture experimental work will focus on the effect of larger 
fields. 
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